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Abstract: NLO corrections to jet cross sections in DIS at HERA are studied, with 
particular emphasis on the two jet final state. High jet transverse momenta are a 
good criterion for the applicability of fixed order perturbation theory. A "natural" 
scale choice is the average k^ of the jets in the Breit frame, which suggest analyzing 
the data in different < k^ > intervalls. 

An important topic to be studied at HERA is the production of multi-jet events in DIS, 
where the expected good event statistics |IJ allows for precision tests of QCD 0. Such tests 
require next-to-leading order (NLO) QCD corrections. Full NLO corrections for one and two- 
jet production cross sections and distributions are now available and implemented in the fully 
differential ep — > n jets event generator MEPJET [[J, which allows to analyze arbitrary jet 
definition schemes and general cuts in terms of parton 4-momenta. A variety of topics can be 
studied with these tools. They include: a) The determination of a s (fiR) from dijet production 
over a range of scales, /zr, b) The measurement of the gluon density in the proton (via jg — > qq), 
c) Associated forward jet production in the low x regime as a signal of BFKL dynamics [|J. 

The effects of NLO corrections and recombination scheme dependences on the 2-jet cross 
section were discussed in Refs. || ||, |6| already for four different jet algorithms (cone, kj>, JADE, 
W) . While these effects are small in the cone and kx schemes, very large corrections can appear 
in the VT-scheme or the modified JADE scheme, which was introduced for DIS in Ref. 0. 

At leading order (LO) the W and the JADE scheme are equivalent. The NLO cross sections 
in the two schemes, however, can differ by almost a factor of two [^|, [5|], depending on the 
recombination scheme and on the definition of the jet resolution mass ((Mfj = (pi + Pj) 2 in 
the W scheme versus Mf, = 2EiEj(l — cos%) defined in the lab frame in the JADE scheme). 
Trefzger and Rosenbauer || find similarly large differences in the experimental jet cross sections 
(which are in good agreement with MEPJET predictions^), when the data are processed with 
exactly the same jet resolution mass and recombination prescription as used in the theoretical 
calculation. The large differences between and within the JADE and W schemes are caused by 
sizable single jet masses (compared to their energy), predominantly for jets in the central part 

-"■The two jet rate in the W scheme (with E recombination) and in the JADE scheme for corrected ZEUS 
data are 18.6 ± 0.7% and 8.6 ± 0.5%, respectively. The corresponding NLO predictions from MEPJET for the 
same kinematics and the same jet definitions are 17.9% and 8.6 % (see T. Trefzger (2|). 



of the detector. Such single jet mass effects first appear in a NLO calculation where a jet may 
be composed of two partons. Clearly, theoretical calculations must be matched to experimental 
definitions and such potentially large single jet mass effects must be taken into account. 

Previous programs || |9| were limited to a.W type algorithm^ and are not flexible enough to 
take into account the effects of single jet masses or differences between recombination schemes. 
In addition, approximations were made to the matrix elements in these programs which are 
not valid in large regions of phase space ||. These problems are reflected in inconsistent 
values for ot a (M|) [ranging, for example, from 0.114 to 0.127 in the HI analysis |1[, (see K. 
Rosenbauer @)], when these programs are used to analyze the data with different recombination 
schemes. Because of these problems, the older programs cannot be used for precision studies 
at NLO in their present form [Kj. In order to reduce theoretical errors, previous analyses 



||T| should be repeated with MEPJET or a similar flexible Monte Carlo program [ IT] . A first 
reanalysis, with MEPJET, of HI data by K. Rosenbauer yields a markedly lower central value, 
a s (M§) = 0.112, which is independent of the recombination scheme (used in both data and 
theory), and the a s (/z|j) extracted from different kinematical bins follows nicely the expectation 
from the renormalization group equation. A similar reanalysis of the ZEUS data has already 
been performed by T. Trefzger, also with MEPJET. 

Single jet mass effects and recombination scheme dependences are fairly small in the cone 
and kx schemes || which, therefore, appear better suited for precision QCD tests. In the 
following, we concentrate on these two and the E recombination scheme. A first issue which 
must be addressed is the dependence of the NLO 2-jet cross section on the renormalization 
scale, fiR, and the factorization scale, \ip. The chosen scale should be characteristic for the 
QCD portion of the process at hand. For dijet invariant masses, rrijj, below Q we are in 
the DIS limit and Q is expected to be the relevant scale. For large dijet invariant masses, 
however, 3> Q, the situation is more like in dijet production at hadron colliders and the jet 
transverse momenta set the physical scale of the process. A variable which interpolates between 
these two limits is the sum of jet kps in the Breit frame 0, J2j ^tU)- Here, {k§(j)) 2 is defined 
by 2E?(1 — cos 9 j P ), where the subscripts j and P denote the jet and proton, respectively. 
J2j &t0) approaches Q in the parton limit and it corresponds to the sum of jet transverse 
momenta, p T , (with respect to the 7*-proton direction) when the photon virtuality becomes 
negligible. We use this "natural" scale for multi-jet production in DIS in the following. 

A good measure of the improvement of a NLO over a LO prediction is provided by the 
residual scale dependence of the cross section. As an example we use the kx algorithm (imple- 
mented in the Breit frame) as described in Ref . [ VI \ . One finds very similar results for the cone 



scheme. Kinematical cuts are imposed on the final state lepton and jets to closely model the HI 
event selection [O. More specifically, we require 10 GeV 2 < Q 2 < 10000 GeV 2 , 0.01 < y < 1, 



0.0001 < x < 1, and an energy cut of E(l') > 10 GeV and a cut on the pseudo-rapidity 
r] = — lntan(6'/2) of the scattered lepton. This r\ cut is Q 2 dependent: —2.794 < r](l') < —1.735 
for Q 2 < 100 GeV 2 and -1.317 < r}(l') < 2.436 for Q 2 > 100 GeV 2 . In addition, we require 
— 1.154 < r](j) < 2.436. The hard scattering scale, E 2 ^, in the kj- algorithm is fixed to 40 GeV 2 
and y cu t = 1 is the resolution parameter for resolving the macro-jets. 

Fig. la shows the scale dependence of the dijet cross section in LO and NLO for the kp 



scheme. The LO (NLO) results are based on the LO (NLO) parton distributions of GRV [|TJ 



2 DISJET H and PROJET || are largely based on the fact that the calculation of the jet resolution mass 
squared, My, can be done in a lorentz invariant way, i.e. as in the W scheme. Only in LO does this agree with 
the JADE definition, defined in the lab frame. 



2000 



1800 



1400 



1300 



"I I I 



a) 



1000 „ 
10" 2 10 




> 

CD 

o 

\ 

XI 

a. 



300 



200 



100 




5 10 15 20 25 30 
<k T (jet)> [GeV] 



10* 



XI 

£: 

<y 10 1 

T3 

\ 
b 



10 u 



10" 




5 10 15 20 25 30 
Q [GeV] 



Figure 1: a) Dependence of the two-jet exclusive cross section in the kr scheme on the scale 
factor £. The solid curves are for fj, R = jj? F = 

€ (Ei ^r(0) 2 > while for the dashed (dotted) 
curves only £ R — £ (£p = £) is varied but £ F = 1/4 (£ R = 1/4) is fixed. Results are shown 
for the LO and NLO calculations, b) NLO < > distribution for the two-jet exclusive cross 
section, c) NLO Q distribution for the four bins in b). 



together with the one-loop (two- loop) formula with five flavors for the strong coupling constant. 
The scale factors £ are defined via 

/4 = & (E kid)? , & = ip (E km 2 • (i) 

i i 

The LO variation by a factor 1.55 is reduced to a 11% variation at NLO when both scales 
are varied simultaneously over the plotted range (solid curves). Also shown is the £ = £_r 
dependence of LO and NLO cross sections at fixed £p = 1/4 (dashed curves) and the £ = £p 
dependence of LO and NLO cross sections at fixed £r = 1/4 (dotted curves). The NLO 
corrections substantially reduce the renormalization and factorization scale dependence. If not 
stated otherwise, we fix the scale factors to £ = £r = £f = 1/4 in the following discussion. 

Let us denote the average k? of the (two) jets in the Breit frame by 

< 14 >= \ ( E >4U))- (2) 

Z 3=1,2 

Fig. lb shows the < k^ > distribution for the NLO 2-jet exclusive cross section in the kx 
scheme. We divide the distribution into four < k^ > bins (suggesting a separate determination 
of a s (< k§, > 2 ) for each). The dependence of the NLO cross section on the scale factor, £, is 
shown in Table 1 for individual bins, and is typically below ±5%. These fairly small theoretical 
uncertainties in the kr algorithm are due to the relatively high value of the hard scattering 
scale, E\ > 40 GeV 2 (or roughly equivalent cuts of p l ^ b ,p1 > 5 GeV on the jets in the cone 
scheme). Thus a precise measurement of a s (< k^ > 2 ) should be possible. 

The Q distributions for the NLO exclusive dijet cross section for these four bins in Fig. lc 
show that even events with very large < k? > are dominated by the small Q 2 region. (The 
dips in the Q distribution around Q = 10 GeV are a consequence of the rapidity cuts on the 



Table 1: NLO (LO) 2-jet exlusive cross sections in pb for the four < > bins and their sum. 
Results are shown for three different choices of the scale factor £ = £r = £p. 

e = i e = 1/4 e = i/ie 



bin 1 
bin 2 
bin 3 
bin 4 



5 GeV < 
10 GeV < 
15 GeV < 
20 GeV < 

sum of bins 



k^p 



< 10 GeV 

< 15 GeV 

< 20 GeV 



881 (821) 
396 (357) 
105 (102) 
63 (68) 
1445 (1348) 



900 (907) 
415 (403) 
107 (118) 
64 (80) 
1486 (1508) 



934 (999) 
433 (461) 
106 (137) 
57 (95) 
1530 (1692) 



scattered lepton, see above). Thus there is a qualitative difference betwen scale choices tied 
to < > versus scales related to Q. One finds that \j? R , /i 2 F = ^Q 2 gives a much larger £ 
dependence for dijet events at NLO than the ones exhibited in Fig. la B. This is the reason 
why scales tied to k? are better suited for QCD analyses of multijet events in DIS. 
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Abstract: NLO corrections to jet cross sections in DIS at HERA are studied, with 
particular emphasis on the two jet final state. High jet transverse momenta are a 
good criterion for the applicability of fixed order perturbation theory. A "natural" 
scale choice is the average kj of the jets in the Breit frame, which suggest analyzing 
the data in different < kj > intervalls. 

An important topic to be studied at HERA is the production of multi-jet events in DIS, 
where the expected good event statistics [I] allows for precision tests of QCD [2]. Such tests 
require next-to-leading order (NLO) QCD corrections. Full NLO corrections for one and two- 
jet production cross sections and distributions are now available and implemented in the fully 
differential ep — >■ n jets event generator MEPJET [3], which allows to analyze arbitrary jet 
definition schemes and general cuts in terms of parton 4-momenta. A variety of topics can be 
studied with these tools. They include: a) The determination of o: s (//r) from dijet production 
over a range of scales, //r, b) The measurement of the gluon density in the proton (via ~fg — > qq) } 
c) Associated forward jet production in the low x regime as a signal of BFKL dynamics [4]. 

The effects of NLO corrections and recombination scheme dependences on the 2-jet cross 
section were discussed in Refs. [3, 5, 6] already for four different jet algorithms (cone, JADE, 
W) . While these effects are small in the cone and kj schemes, very large corrections can appear 
in the VU-scheme or the modified JADE scheme, which was introduced for DIS in Ref. [7]. 

At leading order (LO) the W and the JADE scheme are equivalent. The NLO cross sections 
in the two schemes, however, can differ by almost a factor of two [3, 5], depending on the 
recombination scheme and on the definition of the jet resolution mass ((M^ = (pi + pj) 2 in 
the W scheme versus Mf- = 2EiE 1 (l — cos 6ij) defined in the lab frame in the JADE scheme). 
Trefzger and Rosenbauer [2] find similarly large differences in the experimental jet cross sections 
(which are in good agreement with MEPJET predictions 1 ), when the data are processed with 
exactly the same jet resolution mass and recombination prescription as used in the theoretical 
calculation. The large differences between and within the JADE and W schemes are caused by 
sizable single jet masses (compared to their energy), predominantly for jets in the central part 
of the detector. Such single jet mass effects first appear in a NLO calculation where a jet may 

1 The two jet rate in the W scheme (with E recombination) and in the JADE scheme for corrected ZEUS 
data are 18.6 ± 0.7% and 8.6 ± 0.5%, respectively. The corresponding NLO predictions from MEPJET for the 
same kinematics and the same jet definitions are 17.9% and 8.6 % (see T. Trefzger [2]). 



be composed of two partons. Clearly, theoretical calculations must be matched to experimental 
definitions and such potentially large single jet mass effects must be taken into account. 

Previous programs [8, 9] were limited to a W type algorithm 2 and are not flexible enough to 
take into account the effects of single jet masses or differences between recombination schemes. 
In addition, approximations were made to the matrix elements in these programs which are 
not valid in large regions of phase space [3]. These problems are reflected in inconsistent 
values for a s (Af|) [ranging, for example, from 0.114 to 0.127 in the HI analysis [1], (see K. 
Rosenbauer [2])], when these programs are used to analyze the data with different recombination 
schemes. Because of these problems, the older programs cannot be used for precision studies 
at NLO in their present form [10]. In order to reduce theoretical errors, previous analyses 
[1] should be repeated with MEPJET or a similar flexible Monte Carlo program [11]. A first 
reanalysis, with MEPJET, of HI data by K. Rosenbauer yields a markedly lower central value, 
a s (Mz) = 0.112, which is independent of the recombination scheme (used in both data and 
theory), and the a s (jj 2 R ) extracted from different kinematical bins follows nicely the expectation 
from the renormalization group equation. A similar reanalysis of the ZEUS data has already 
been performed by T. Trefzger, also with MEPJET. 

Single jet mass effects and recombination scheme dependences are fairly small in the cone 
and kj schemes [3] which, therefore, appear better suited for precision QCD tests. In the 
following, we concentrate on these two and the E recombination scheme. A first issue which 
must be addressed is the dependence of the NLO 2-jet cross section on the renormalization 
scale, //r, and the factorization scale, fip. The chosen scale should be characteristic for the 
QCD portion of the process at hand. For dijet invariant masses, m rjl below Q we are in 
the DIS limit and Q is expected to be the relevant scale. For large dijet invariant masses, 
however, m rj ^> Q } the situation is more like in dijet production at hadron colliders and the jet 
transverse momenta set the physical scale of the process. A variable which interpolates between 
these two limits is the sum of jet kjs in the Breit frame [6], J2 3 &x(i)- Here, (kj(j)) 2 is defined 
by 2 E 2 (l — cos 6jp), where the subscripts j and P denote the jet and proton, respectively. 
J2j kxij) approaches Q in the parton limit and it corresponds to the sum of jet transverse 
momenta, pj,, (with respect to the 7*-proton direction) when the photon virtuality becomes 
negligible. We use this "natural" scale for multi-jet production in DIS in the following. 

A good measure of the improvement of a NLO over a LO prediction is provided by the 
residual scale dependence of the cross section. As an example we use the kj algorithm (imple- 
mented in the Breit frame) as described in Ref. [12]. One finds very similar results for the cone 
scheme. Kinematical cuts are imposed on the final state lepton and jets to closely model the HI 
event selection [13]. More specifically, we require 10 GeV 2 < Q 2 < 10000 GeV 2 , 0.01 < y < 1, 
0.0001 < x < 1, and an energy cut of E(l') > 10 GeV and a cut on the pseudo-rapidity 
n = — lntan(#/2) of the scattered lepton. This i] cut is Q 2 dependent: —2.794 < f](l') < —1.735 
for Q 2 < 100 GeV 2 and -1.317 < ??(/') < 2.436 for Q 2 > 100 GeV 2 . In addition, we require 
— 1.154 < r)(j) < 2.436. The hard scattering scale, E^ } in the kj algorithm is fixed to 40 GeV 2 
and y cut = 1 is the resolution parameter for resolving the macro-jets. 

Fig. la shows the scale dependence of the dijet cross section in LO and NLO for the kj 
scheme. The LO (NLO) results are based on the LO (NLO) parton distributions of GRV [14] 

2 DISJET [9] and PROJET [8] are largely based on the fact that the calculation of the jet resolution mass 
squared, Mfj, can be done in a lorentz invariant way, i.e. as in the W scheme. Only in LO does this agree with 
the JADE definition, defined in the lab frame. 
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Figure 1: a) Dependence of the two-jet exclusive cross section in the kj scheme on the scale 
factor ^. The solid curves are for n\ = [i 2 F = £ kj(i)) 2 , while for the dashed (dotted) 

curves only £r = £ (£p = £) is varied but £p = 1/4 (£r = 1/4) is hxed. Results are shown 
for the LO and NLO calculations, b) NLO < kj > distribution for the two-jet exclusive cross 
section, c) NLO Q distribution for the four bins in b). 



together with the one-loop (two-loop) formula with five flavors for the strong coupling constant. 
The scale factors £ are defined via 

/4 = ME^x(0) 2 , /4 = ME^(0) 2 - (i) 

i i 

The LO variation by a factor 1.55 is reduced to a 11% variation at NLO when both scales 
are varied simultaneously over the plotted range (solid curves). Also shown is the £ = £r 
dependence of LO and NLO cross sections at hxed £p = 1/4 (dashed curves) and the £ = £p 
dependence of LO and NLO cross sections at hxed £r = 1/4 (dotted curves). The NLO 
corrections substantially reduce the renormalization and factorization scale dependence. If not 
stated otherwise, we hx the scale factors to £ = £r = £p = 1/4 in the following discussion. 

Let us denote the average kj of the (two) jets in the Breit frame by 

< 14 >= \ ( E *r0'))- (2) 

Z j = l,2 

Fig. lb shows the < kj > distribution for the NLO 2-jet exclusive cross section in the kj 
scheme. We divide the distribution into four < kj > bins (suggesting a separate determination 
of a s (< kj > 2 ) for each). The dependence of the NLO cross section on the scale factor, £, is 
shown in Table 1 for individual bins, and is typically below ±5%. These fairly small theoretical 
uncertainties in the kj algorithm are due to the relatively high value of the hard scattering 
scale, Ej > 40 GeV 2 (or roughly equivalent cuts of p l j b ,Pj > 5 GeV on the jets in the cone 
scheme). Thus a precise measurement of a s (< kj > 2 ) should be possible. 

The Q distributions for the NLO exclusive dijet cross section for these four bins in Fig. lc 
show that even events with very large < kj > are dominated by the small Q 2 region. (The 
dips in the Q distribution around Q = 10 GeV are a consequence of the rapidity cuts on the 



Table 1: NLO (LO) 2-jet exlusive cross sections in pb for the four < k T > bins and their sum. 
Results are shown for three different choices of the scale factor £ = £r = £p. 





£ = 1 


€ = 1/4 


€ = 1/16 


bin f: 


5 GeV < < k$ 


> 


< fO GeV 


881 (821) 


900 (907) 


934 (999) 


bin 2: 


fO GeV < < k$ 


> 


< 15 GeV 


396 (357) 


415 (403) 


433 (461) 


bin 3: 


15 GeV < < k$ 


> 


< 20 GeV 


105 (102) 


107 (118) 


106 (137) 


bin 4: 


20 GeV < < k T 


> 




63 (68) 


64 (80) 


57 (95) 




sum of bins 






1445 (1348) 


1486 (1508) 


1530 (1692) 



scattered lepton, see above). Thus there is a qualitative difference betwen scale choices tied 
to < k T > versus scales related to Q. One finds that //^, fj, F = £Q 2 gives a much larger £ 
dependence for dijet events at NLO than the ones exhibited in Fig. la [3]. This is the reason 
why scales tied to k F are better suited for QCD analyses of multijet events in DIS. 
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